Direct absorption imaging of ultracold polar molecules 
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We demonstrate a scheme for direct absorption imaging of an ultracold ground-state polar molec- 
ular gas near quantum degeneracy. A challenge in imaging molecules is the lack of closed optical 
cycling transitions. Our technique relies on photon shot-noise limited absorption imaging on a 
strong bound-bound molecular transition. We present a systematic characterization of this imaging 
technique. Using this technique combined with time-of-flight (TOF) expansion, we demonstrate the 
capability to determine momentum and spatial distributions for the molecular gas. We anticipate 
that this imaging technique will be a powerful tool for studying molecular quantum gases. 

PACS numbers: 67.85.-d, 37.10.Pq, 33.20.-t 



The field of ultracold polar molecules has gained a 
forefront interest in recent years In contrast to the 
contact interaction between atoms, the electric dipolar 
interactions between polar molecules are long range and 
anisotropic. Together with their complex internal struc- 
ture, polar molecules are promising candidates for in- 
vestigating strongly correlated many-body systems @, 
[!§], quantum information ultracold quantum chem- 
istry [j| @] , and precision measurements J?} ■ The recent 
breakthrough in the production of a near quantum de- 
generate gas of KRb molecules [|[ has paved the way 
for studying ultracold polar molecules with high phase- 
space density. Chemical reactions controlled by quantum 
statistics of the fermionic KRb molecules have recently 
been observed Q. Dipolar interaction effects on these 
molecular collisions have also been studied 0. 

Convenient and powerful detection methods are at the 
heart of successful research with ultracold atoms flfij . 
Optical cycling transitions make it possible for each atom 
to scatter many photons, which greatly enhances the 
signal-to-noise ratio (SNR) for optical absorption imag- 
ing. Images of the cloud can be taken with or without 
TOF expansion and these can be used to obtain the spa- 
tial and momentum distributions of the trapped gas. In 
the case of molecules, the wavef unction overlap between 
ground and excited states (Franck-Condon factor (FCF)) 
is typically much less than unity. The absence of cycling 
transitions then directly implies that the maximum num- 
ber of photons scattered from the molecular target basi- 
cally equals to the number of molecules. Hence, detec- 
tion of a small number of cold molecules formed, for ex- 
ample, via photo-association typicall y re lies on very sen- 
sitive ionization detection methods [ll|. Laser-induced 
fluorescence detection on decelerated molecules has also 
been successfully employed [12, ll2| . Although these tech- 
niques are powerful, it is challenging to make their com- 
plex setups compatible with the ultrahigh vacuum sys- 
tems required for quantum gas experiments. In addi- 
tion, these methods need pulsed lasers to excite or ion- 
ize the molecules, and velocity imaging techniques would 



need to be implemented (14| to obtain momentum infor- 
mation. In our previous work @, 0] , ground-state KRb 
molecules were detected using a complex approach of co- 
herent transfer of the molecules to a weakly bound state 
near a Fano-Feshbach resonance followed by atomic ab- 
sorption imaging. Because of the tiny binding energy, 
these molecules could be photodissociated and detected 
using light on the atomic cycling transitions HI, Hi ]. 

In this Letter, we demonstrate imaging of a high phase- 
space-density sample of ground-state polar molecules us- 
ing direct optical absorption. The ima ging is accom- 



plished via an open optical transition 13, LL8[ . The opti- 
cal density is a dynamic quantity here since it decreases 
rapidly during the probe time as the population decays to 
many dark states. In addition, the absorption cross sec- 
tion for the molecule is usually at least one order of mag- 
nitude smaller than those for atoms, due to the low FCF 
and rotational branching. Here we implement absorption 
imaging on one of the strongest molecular bound-bound 
transitions for KRb. By counting the number of absorbed 
photons via a CCD camera, we obtain molecular imaging 
with a high SNR. The highest SNR is achieved by using a 
minimum number of photons for imaging the molecules 
with photon shot-noise limited detection. The demon- 
strated technique can be used to determine momentum 
and spatial distributions for the molecular gas and al- 
lows for flexible molecular imaging at arbitrary external 
electric (E) and magnetic (B) fields. 

Our experimental procedure to produce ultracold 
40 K 87 Rb molecules in the lowest vibrational (v" — 0) and 
rotational (N" — 0) level of the X 1 E + electronic ground 
state is described in Ref. @- Near quantum-degenerate 
40 K and 87 Rb atoms are prepared in a pancake-shaped 
crossed-beam optical dipole trap with an aspect ratio 
of 6 between the tighter confined vertical direction z 
and the horizontal direction x. After creating weakly 
bound 40 K 87 Rb Feshbach molecules in a high vibrational 
state, we coherently transfer the molecular population to 
a single hyperfine state of the ro-vibronic ground-state 
manifold via stimulated Raman adiabatic passage (STI- 
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FIG. 1: (Color online) (a) The ab initio electronic potentials 
for KRb, showing the direct absorption detection scheme for 
X 1 ^ molecules in the v" = 0, N" = level. The probe light 
drives a transition to the v' = 3 level of the l 1 !! state. Inset: 
Schematic diagram of the open two-level model described in 
the text, (b) An example population depletion lineshape for 
this transition taken with 10 ps probe pulse duration, (c) 
Decay of the ground-state population as a function of the 
probe pulse duration. The probe intensity is ~0.7 mW/cm 2 . 
The solid curves are fits based on the open two-level model. 



RAP) [8j, [l9(. Our standard method to detect these 
ground-state molecules has been to transfer them back 
into the weakly bound Feshbach state, followed by detec- 
tion with the usual atomic absorption imaging method. 
Because of the coherent transfer process, this Feshbach 
imaging scheme can only detect molecules in a single hy- 
perhne state. It is also very sensitive to external helds. 

To image the ground-state molecules directly, we use 
a transition at A = 658.4 nm, from the absolute ground 
state to the vibrational v 1 = 3, rotational J' = 1 level 
of the l 1 !! state, as shown in Fig. [T] (a). The l 1 !! state 
has been studied in_detail for a different isotope combi- 
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nation, cU K f "Rb [20| . and the selected transition is one 
of the strongest in the KRb molecule 2l|. Using the 
available potential (2fl|, we calculate the relevant rovi- 
brational levels for 40 K 87 Rb. The measured transition 
frequency is within 30 GHz of the prediction. Shown in 
Fig.Q](b) is a typical molecular resonance lineshape taken 
at a relatively low intensity of ~0.7 mW/cm 2 . The x-axis 
is the probe laser frequency minus the measured transi- 
tion frequency of 455219 GHz. The number of remain- 
ing molecules is measured using the Feshbach imaging 
method described above. The nature of an open transi- 
tion can be seen in Fig. [T] (c), where we show the loss of 
ground-state molecules as a function of the probe time. 
For these data, the probe laser frequency is on resonance. 

To calibrate the transition strength, we simultaneously 
fit several lineshapes and decay curves obtained for a 
range of intensities to an open two-level model, which is 
shown schematically in the inset of Fig. 1 (a). This model 
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FIG. 2: (Color online) (a) An absorption image (228 frni x 
171 lira) of 39,000 ground-state molecules after 2 ms TOF. 
The false color indicates CCD counts, taken for the best SNR 
with ~75% of the molecules detected, (b) Peak SNR as a 
function of n p hoton- Filled circles: experimentally measured 
peak SNR; Solid curve: fit to the model described in the text. 
The measured A r mo i ecuie = 950 per pixel at the cloud center. 



is described by the optical Bloch equations (OBEs): 

i 

Pgg = Pee + ^(Peg ~ Pge) 

1 

Pee = -IPee + ^(Pge ~ Peg) 
Pge = -i^+i^Pge^ ^^(Pee- Pgg)- (1) 

Here, the density matrix elements p gg and p ee are the 
fractional population in the ground and excited states, 
respectively. The other two elements p eg — p* e are re- 
lated to the coherence. The laser frequency detuning 
from the transition is A. The probe Rabi frequency is 
given by f2 = dS /H, where d is the electronic transition 
dipole moment including the FCF, S is the probe field 
and H is the Planck's constant over 2ir. As the calcu- 
lated FCF is 0.21 [2(| and the rotational branching ratio 
is 2/3 according to the Honl-London coefficients for the 
probe transition, the branching ratio r\ is 0.14. Initially, 
all molecules are prepared in the ground state, p gg = 1 
and p ee = 0. At t > 0, p gg + p ee < 1 due to the de- 
cay rate (1 — 77)7 from the excited state to dark states. 
The excited population also decays back to the ground 
state with a smaller rate 777. Shown by the solid curves in 
Fig.QJb) and (c) are the fits of the data to the OBEs with 
7 and n as free parameters. From the fits, we determine 
7 = 2?r x 22(7) MHz and ft = 2tt x 0.8(1) MHz. 

A typical direct absorption image of ground-state 
molecules is shown in Fig. [2] (a). This image is ob- 
tained using a CCD camera. In one experimental cycle, 
two images are taken, one with molecules and the other 
with just the probe light for reference. By subtracting 
these two images, we obtain the molecular absorption 
image. Furthermore, the subtraction eliminates most of 
the technical background noise, which arises from spa- 
tial variations in the probe beam intensity. Taking into 
account the CCD calibration (including a quantum ef- 
ficiency of 0.9 electron/photon and a digitizing gain of 
1.1 electron/count), an 11% loss from the imaging opti- 



cal elements, and the braching ratio r\ = 0.14, we need 
to multiply by a factor of 1.2 to convert CCD counts 
to molecule number. For each image, we perform a sur- 
face Gaussian fit to obtain the peak signal, N pea f., in 
CCD counts per pixel and sizes <r Xl a z in number of pix- 
els. The molecule number can then be obtained from 
2TTN pea k(T x a z x 1.2. The detected number of molecules 
by direct imaging is typically 15% lower than that mea- 
sured by the Feshbach imaging method. To extract the 
ground-state molecule number from the Feshbach imag- 
ing method, we have to consider the STIRAP efficiency, 
the Feshbach molecule detection efficiency, and optical 
density saturation. Thus, the number estimate from di- 
rect imaging has smaller uncertainties. In addition, in 
contrast to conventional absorption imaging, the open- 
transition imaging eliminates the need for the exact value 
of the absorption cross section. 

Since the total signal size is limited to the molecule 
number, it is crucial to understand the detection SNR 
and optimize it. We have verified that the SNR in our 
experiment is limited by shot noise, which scales with the 
square root of the photon number. Under the condition 
of weak probing (SI << 7), the fractional ground-state 
population follows p gg = e~ n t ( 1_, ')/T. Here t is the probe 
pulse duration and fi 2 t = 47r j-f-y^n p hotom with 

Nphoton 

the average photon number per CCD pixel at the molec- 
ular cloud position. The area per pixel A is 5.7 /xm x 
5.7 /zm at the cloud position and eo is the vacuum per- 
mittivity. As each molecule can scatter 1/(1 — rf) pho- 
tons before decaying into a dark state, we express the 
signal size (total number of scattered photons) per pixel, 
N signal, as a function of n photon , 

N signal = N ™^le {l _ e ^^l nphotonl (2) 

where N mo i ecu i e is the initial number of molecules per 
pixel. The photon shot-noise-limited SNR is then 

SNR = 0.8N slgnal /^/2 x 0.8 x n photon , (3) 

where the factors of 2 in the denominator arises from the 
fact that we subtract two frames to obtain one absorp- 
tion image. The net detection efficiency is 0.8 (including 
losses from optics and quantum efficiency of the CCD). 

Fig. [5] (b) shows the peak SNR as a function of n p hoton 
where n p h oton is varied by using different t at a constant 
probing intensity. The peak SNR is determined from the 
absorption image by comparing the measured number of 
counts per pixel at the cloud center to the standard de- 
viation of the measured number of counts per pixel far 
from the cloud center. The solid line in Fig. [5] (b) shows 
a fit to the model described by Eqns. ([2]) and ([3]) with d 
the only free parameter. The initial number of molecules 
per pixel at the cloud center, N mo iecuie = 950, is deter- 
mined from the data. The fit yields d = 3.5(6) Debye, 
consistent with measurement of ref. [20]. The SNR ini- 
tially increases as n p hoton is increased; however, for large 
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FIG. 3: (Color online) (a) Measurement of the molecular gas 
temperature by TOF expansion followed by direct absorption 
imaging. The trap frequency along the x and z direction 
is 32 Hz and 195 Hz, respectively. Filled circles and open 
triangles are the rms cloud widths in the x and z directions, 
respectively. Solid curves are fits to a ballistic expansion, (b) 
Measurement of the molecular two-body loss rate /3, for an 
initial temperature of 320 nK, using the direct imaging (open 
circles) and Feshbach imaging (filled squares) methods. Solid 
curves are fits to a two-body loss model. The extracted values 
of ft agree within the experimental uncertainty. 

n P hotom the growth of photon shot noise dominates the 
signal increase. From Eqns. ^ and ([3]), the best SNR 
occurs when 71% of the population is depleted by the 
probing light. In a typical experiment, we use n p hoton 
slightly greater than that for the optimal SNR to ensure 
that all molecules are depleted. This allows us to accu- 
rately count the total number of molecules. 

Fig. |3] (a) demonstrates the use of direct absorption 
imaging to measure the temperature of a molecular cloud 
using TOF expansion following sudden release from the 
optical trap. We measure molecular cloud sizes for ex- 
pansion times up to 9 ms. As the molecular cloud ex- 
pands, the detection SNR decreases. For longer expan- 
sion times, we average several successive images to im- 
prove the SNR. From the measurements in Fig. [3] (a) , we 
extract a mean cloud temperature of 220(13) nK, which 
corresponds to ~1.4 times the gas Fermi temperature. 

We have also taken a lifetime measurement of the 
trapped molecules to study inelastic collisions. In 
Ref. @, the observed molecular two-body loss, detected 
by Feshbach imaging, was due to chemical reaction of 
KRb + KRb. The fact that the molecules were pre- 
pared in the lowest energy state, including the hyper- 
fine structure [h| , led to the unambiguous observation of 
chemical reactions. If the KRb molecules are instead pre- 
pared in an excited hyperfine state, inelastic hyperfine- 
state-changing collisions are also a possible loss mecha- 
nism. With the Feshbach imaging method, we cannot 
distinguish between the two possible loss channels be- 
cause the imaging is only sensitive to a single hyperfine 
state. In the direct imaging method, the probe transition 
linewidth is much larger than the entire span of ~ 3 MHz 
of the ground-state hyperfine manifold at 550 Gauss. The 
direct imaging thus provides a simultaneous detection of 
all the molecules regardless of their distribution among 
different hyperfine states. Hence, a comparison of direct 
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FIG. 4: (Color online) Measurement of the Zeeman (a, black 
filled circles) and Stark (b, blue filled circles) shifts of the 
optical transition used for direct absorption imaging. Open 
triangles correspond to excited state Stark shifts obtained af- 
ter subtracting the known Stark shift of the ground state []| ; 
Solid curves are quadratic fits. 

imaging and Feshbach imaging could allow us to check 
if nuclear spin flips to other hyperfme states occur in an 
ultracold molecule gas. Fig. [3] (b) presents exactly such 
a comparison. The data taken here are for ground-state 
molecules prepared in an excited hyperfme state. At a 
temperature of 320 nK, the loss rate coefficients obtained 
using the two imaging techniques are in good agreement. 
While this agreement is consistent with there being no 
appreciable rate for hyperfine-state changing collisions, 
at the large B field used here, it is also possible that 
spin-flipped molecules might acquire a sufficient kinetic 
energy to leave the trap before imaging. 

In future work, it will be interesting to explore 
molecule loss rates for different B or E fields. Here, 
the direct molecular imaging method is advantageous be- 
cause it can operate under any combinations of E and 
B fields. Besides investigating the possibility of nuclear 
spin flips, the technique could also be useful in looking for 
molecular scattering resonances induced by E [22| or B 
fields. In addition, an important feature of direct imag- 
ing is the capability of imaging the molecular gas in situ 
under the presence of strong dipolar interactions. 

In Fig. 21 we show that the probing transition for di- 
rect molecule imaging is relatively insensitive to B and 
E fields. The v' = 3, J' = 1 level of the ^ state has 
first-order Zeeman and Stark shifts due to a half Bohr 
magneton magnetic dipole moment and a non-zero elec- 
tric dipole moment. To minimize variation in the reso- 
nant frequency for probing over a broad range of E and 
B fields, we use the Mj = sublevel as the excited state 
for detection. In doing so we eliminate the first-order 
Zeeman and Stark shifts in the excited state as they are 
both proportional to Mj. Figure 0] shows the observa- 
tion of relatively small (a) Zeeman and (b) Stark shifts 
over a large range of B and E fields. The quadratic 
dependences of the small shifts are indicative of the sec- 
ond order Zeeman and Stark effects, respectively. The 
ground state has a negligible magnetic dipole mo- 
ment, due only to the nuclear spin. Thus the observed 



Zeeman shift in Fig. [4] (a) arises solely from the excited 
state. The Stark shift of the ground state was measured 
previously Q . The excited-state Stark shift can be deter- 
mined by subtracting the ground-state Stark effect from 
the total observed shift (Fig. g] (b)). 

In conclusion, we have demonstrated direct absorption 
imaging of an ultracold molecular gas. Whereas the pre- 
vious technique of coherent transfer followed by absorp- 
tion imaging of Feshbach molecules is highly state selec- 
tive, direct imaging allows us to probe molecules indepen- 
dent of their hyperfme state. The detection scheme also 
eliminates the extreme sensitivity on B and E fields, en- 
abling investigation of dipolar effects under various com- 
binations of external fields. The detection scheme should 
be universal for all alkali dimers due to their spectro- 
scopic similarities. We note that a non-destructive phase- 
contrast imaging technique [23| might be particularly ad- 
vantageous for open structured molecules. 
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